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ABSTRACT 

We investigate the luminosity and cooling of highly magnetised white dwarfs. We 
consider white dwarfs with electron-degenerate core and nondegenerate surface layers 
where cooling occurs by diffusion of photons. We find the temperature and density 
profiles in the surface layers or envelope of white dwarfs for radially constant and 
varying magnetic fields by solving the magnetostatic equilibrium and photon diffusion 
equations in a Newtonian framework. We also obtain the properties of white dwarfs 
at the core-envelope interface, when the core is assumed to be practically isothermal 
due to large thermal conductivity. With the increase in magnetic field, the interface 
temperature and density are found to be increasing. While the interface radius also 
increases with the increase in magnetic field when the field is hypothesised to be con¬ 
stant throughout the star, the interface radius decreases for varying fields. However, 
for white dwarfs having fixed interface radius or interface temperature, we find that the 
luminosity significantly decreases, falling in the range ~ 10“® — 10“^^ solar luminosity, 
with the increase in magnetic field strength at the interface and hence envelope, in the 
corresponding range 10® —10^^ G, in particular for the varying magnetic fields which 
are expected to be more realistic. This is remarkable as it argues for magnetised white 
dwarfs to be dimmer and be practically hidden in the H-R diagram. We also find the 
cooling rates corresponding to these luminosities. Interestingly, the decrease in tem¬ 
perature with time, for the fields under consideration, is not found to be appreciable 
— at most by a factor of two and that is also for the constant field cases. 

Key words: stars: magnetic fields - white dwarfs - stars: luminosity function, mass 
function - MHD - stars: massive 


1 INTRODUCTION 


One of the most puzzling observations in high energy astrophysics in the last decade o r so is that of the over-luminous type l a 
supernova. More than a dozen such supernovae have been observed since 2006 (e.g. Howell et al. 20061 : Scaizo et al. 2O10l '). 
whose significantly high luminosities can only be justified by invoking very massive progenitors, of mass M 2Mq where Mq 
is solar mass. However, the Chandrasekhar-limit prohibits the existence of such a massive progenitor white dwarf. Proposed 
models explaining the se highly super-Chandrasekhar progenitors include rapidly (and differen tially) rotating white dwarfs 


( Yoon fc Langer 200J) and binary evolution of accreting differentially rotating white dwarfs ( Hachisr] lOSfj ). Another set 


of proposals which has recently brought the issue of super-Chandrasekhar white dwarfs into the limelight relates to highly 
magnetised white dwarfs. In a series of papers, the main message of this work has been that the enormous power of a 
magnetic field — irrespective of it s natu r e of origin: quantum (due to const ant super-strong field: iKundu fc Mukhopadhvav 
2OI2I : Das fc Mukhopadhvav 2012I . 20131 : Das. Mukhopadhvav fc Rao 2013h . classical and/or general relativistic (due to a 
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varying strong field exerting magnetic pressure and tension: IPas fc Mukhopadhyay 2014a . 2015 1 — is capable of revealing 
significantly super-Chandrasekhar white dwarfs. Remarkably, unlike other proposals, this work also adequately predicts the 
required mass-range 2.1 — 2.8 Mq of the progenitors in order to explain the set of over-luminous type la supernovae. 

The impact of high magnetic fields not only lies in increasing the limiting mass of white dwarfs but it is also expected 


to change other properties including luminosity, temperature, cooling rate etc. Indeed, observations llValvavin et al.l l2014ll 


argue for the suppression of cooling in white dwarfs due to strong magnetic fields. This is in particular due to suppression of 
convection over the entire surface of magnetised white dwarfs, which in turn hinders their cooling evolution. Therefore, highly 
magnetised white dwarfs may have lower surface temp erature, and hence, luminos i ty. Nevertheless, poloidally dominate d 
magnetised white dwarfs are shown to be smaller in size ( Das fc Mukhopadhvavl[2015l : Subramanian fc Mukhopadhvavll2015 l. 

which can account for their lower luminosity. _ 

Although magnetised white dwarfs were attempted to explore earlier (e.g. Ostriker fc Hartwick 19681 : Adam 


author concentrated on the effects of magnetic fields on the internal properties of white dwarfs, e.g. thermal conduction, cooling 
rate etc. as well as luminosity, even thou gh sometime s the chosen field was well above the critical field. Be = 4.414 x 10^^ 
G for Landau effects to arise (see, e.g., Adam 1986h . The same is true for the exploration of magnetised white dwarfs 
by the Mukhopadhyay-group, mentioned above. Such an investigation has been more important for the magnetised white 
dwarfs proposed by the Mukhopadhyay group, because the underlying central fields could be stronger, sometimes > 5 x 10^® 
G. Nevertheless, super-Ghandrasekhar, magnetised white dwarfs were also explored with relatively weaker central fields, e.g. 
~ 5 X 10^"^ G, which do not reveal any quantum mechanical effects, but the underlying magnetic pressure gradient, determined 
by the field geometries and profiles, is respo nsible for making the mass super-Chandrasekhar (IDas fc Mukhopadhvavll2014al . 


2 OI 5 I : Subramanian fc Mukhopadhyay 2015 1. 


In the present paper, we study the cooling of magnetised white dwarfs and calculate the corresponding lumi nosities . While 
the cooling of white dwarfs is not a completely resolved issue, it has been investigated since 1950s, e.g., when Mestel ( 1952ll 
attempted to understand th e source of energy of white dwarfs and to estimate the ages of observed white dwarfs. Subsequently, 


Mestel fc RudermanI (119671 1 also explored the cooling of white dwarfs and found them to be radiating at the expense of their 


thermal energy. They also showed, by including the Goulomb contribution to the thermal energy, that the lifetimes of white 
dwarfs are twice as large as estimated earlier. The evolution and cooling of low mass white dwarfs, beginning as a bright 
central star to the stage of crystallization after about 10 billion years, were also addr essed by Tulmkov fc YungelsonI ( IQQOl l 
who argued that the similarity o f a modern coo ling curve then to the one predicted bv iMestell ( 1952l l is the consequence of a 
series of accidents. Furthermore, Hansen ( 1999l l reviewed the physics of cool white dwarfs, in particular their impact in order 
to extract valuable information about the early history of our galaxy. 

The above works either did not consider the effects of magnetic field or the fields embedding the star were assumed to 
be too weak to have any practical effects. On the other hand, the field of magnetised white dwarfs considered by our group 
(and some others), as mentioned above, is high enough compared to that of all previous works which addressed the cooling of 
white dwarfs. Hence, here we plan to explore the cooling and luminosity of magnetised white dwarfs proposed by us earlier. 

This paper is organized as follows. In §2, we recapitulate how the density and temperature profiles at the interface can be 
obtained for non-magnetised white dwarfs, in addition to computing their cooling timescales. In §3, we include the contribution 
of the magnetic field to the pressure and density of white dwarfs and compute the resultant density and temperature profiles 
for different luminosities and magnetic field strengths. Subsequently, in §4, we consider white dwarfs having either fixed 
interface radius or fixed interface temperature and evaluate their luminosities for increasing field strength. In §5, we compute 
the cooling rates of magnetised white dwarfs for the cases discussed in §4. Next, we discuss the implications of our results for 
non-magnetised and magnetised white dwarfs in §6. Finally, we conclude with a summary in §7. 


2 TEMPERATURE PROFILE AND COOLING FOR A NON-MAGNETISED WHITE DWARF 


In this section, we recapitulate how the luminosity-temperature relation and the temperature profile can be obtained using 
the equations of hydrostatic equilibrium and photon diffusion for a non-magnetised non-relativistic white dwarf having a 
nondegenerate ideal gas envelope. Once we know the luminosity and the specific heat in terms of tem perature, we can find 
the cooling rate of the white dwarf. We use the same notation as used by Shapiro fc Teukolskv ( IQSfil l. 


2.1 Density and temperature profiles in the surface layers 

The electrons in the interior of a white dwarf are completely degenerate and have a large mean free path. This results in a 
high thermal conductivity and a uniform temperature in the inner regions (core) of the white dwarf. Nondegenerate surface 
layers, which form the envelope and are in “radiative equilibrium” with each other and cover the isothermal interior of the 
white dwarf. By radiative equilibrium, we mean that the matter is in thermodynamic equilibrium locally and energy flux is 
carried outward by photon diffusion. 
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Figure 1. Variation of density with temperature for different 
luminosities: 1 O~^L 0 (dashed line), 1 O“‘^L 0 (dotted line) and 
1 O“^L 0 (dot-dashed line), when R = 5000 km and Tg = 5000 
K. 



Figure 2. Variation of radius with temperature for luminosity 
L = 1 O“®-L 0 , when R = 5000 km and Ts = 5000 K. 


Using the equations for photon diffusion and hydrostatic equilibrium, we obtain, 

dP 4ac 47vGm{r) 

It ^ 




3 kL ^ 

In the above equation, P is the pressure, T is the temperature, a is the radiation constant, c is the speed of light in vacuum, 
G is Newton’s gravitational constant, m(r) is the mass enclosed within radius r, n is the opacity of the stellar material and 
L is the luminosity. The opacity can be approximated by Kramer’s formula 

K = KopT~^'^, (2) 

where kq = 4.34 x + X) cm^g“^ and p, X and Z are the density of the stellar material, mass fractions of hyd rogen 

and heavy elements (elements other than hydrogen and helium) in the stellar interior, respectively I Schwarzschildl kssl l . The 
opacity is due to the bound-free and free-free transitions of electrons I Shapiro fc Teukolskv 1983l l. 

Surface layers are expected to be much thinner as compared to the radius of the white dwarf and therefore, m{r) « M 
where, M is the mass of the white dwarf. Henceforth, the solar luminosity will be denoted as Lq and the terms “surface layers” 
and “envelope” will be used interchangeably. In Eqn. (1), we substitute p using the equation of state (EoS) of nondegenerate 
matter (ideal gas) and integrate with the boundary condition P = 0 at T = 0 to obtain 

f 2 iac A-kGM pm^ \ ^ 3.25 , . 

^ \8.5 3 koL kB ) ’ ^ ’ 

where is the atomic mass unit, p is the average mass of particles in units of and fcs is the Boltzmann constant. The 
above relation between p and T is valid only in the surface layers and breaks down for the degenerate core of the white dwarf. 
At the interface between the degenerate core and the nondegenerate envelope, the variation of density with temperature can 
be obtained, by equating the electron pressure on both sides, as 

p* « (2.4 X 10"® g cm-®) (4) 

where pe is the mean molecular weight per electron in units of m^. For a typical white dwarf, X = 0, the mass fraction of 
helium Y = 0.9, Z = 0.1 and M ~ Mq. Using pe « 2, p « 1.4 and combining Eqns. (3) and (4), the luminosity at the 
interface can be written in terms of temperature as 

M 


L«(2x 10® ergs-®)—Tl 
Mg) 


( 5 ) 


Hence, the interface temperature of a white dwarf can be determined if L, M and the composition are known. For observed 
white dwarfs, L « 10-® — 10~^Lq which corresponds to T* « 10® — 10^ K and p* « 10^ — 10® gem-®. Substituting the 


expression for p in the photon diffusion equation and integrating across the envelope, we obtain 


T,-Ts = 


1 prriij, GM 


- 1 


( 6 ) 


4.25 ks R 

where Ts is the temperature at the surface, with r* and R being the radii of the interface and the surface, respectively. 

We conside r white dwarfs with M = Mq. Using the Chandrasekhar’s relation for white dwarfs, this corresponds to radius, 
R = 5000 km ( Chandrasekhar Il931all^ . We will consider this set of mass and radius throughout the paper. However, the 
results presented here practically do not change for other values of R (e.g., in the range 500 — 5000 km) and M, unless Ts 
is as high as 10® K which we do not observe. Figures 1 and 2 show the variations of density and radius, respectively, with 
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temperature in the non-degenerate envelope. Figure 1 shows that the density at a given temperature, and hence radius, is 
suppressed with increasing luminosity. From Eqn. (6), it can be seen that the T — r relation is independent of luminosity. 

Hence, we obtain the relation to be of straight lines with the same slope for different values of luminosity, as shown in 
Fig. 2 for L = 10“®I/q. As T* is different for different luminosities, the corresponding T — r lines originate from different 
temperatures at the interface. 

With the increase of L, T* increases whereas r* decreases. For T* « 10® — 10^ K and Ts « 10® — 10^ K, 

^ 5^ 0.99. (7) 

Corresponding to 10® ^ T* ^ 10^ K, we have p* Si 1500 gcm“® <C Pc, where pc is the central density of the white dwarf. 
Table 1 shows the variations of T*, p, and r, as L changes for a given value of Ts and R. We see that as L increases, T, 
increases whereas r* decreases. 

Table 1: T,, p* and r, for different L, when Ts = 5000 K and R = 5000 km 


L (inLo) 

r. (inK) 

p* (ingcm“®) 

r* (in R) 

1 Ar 1 r*—R ^ > 

10 "® 

2.332 X 10® 

170.722 

0.9978 

2.1097 

5 X 10”® 

3.69346 X 10® 

340.289 

0.9965 

2.11616 

10 "^ 

4.50238 X 10® 

457.995 

0.9958 

2.11841 

5 X 10~'‘ 

7.13096 X 10® 

912.892 

0.9933 

2.12588 

10 ”® 

8.69273 X 10® 

1228.66 

0.9918 

2.12935 

5 X 10"® 

1.37677 X 10^ 

2449.01 

0.9871 

2.14039 

10 "^ 

1.6783 X 10^ 

3296.13 

0.9844 

2.14662 


Hence, as the luminosity of a non-magnetised white dwarf increases, the interface shifts inwards and the degenerate region 
shrinks in volume. However, for the observed range of luminosities, the decrease in volume of the degenerate region is quite 
small. Also |Ar/Ar| increases with luminosity but does not vary appreciably with luminosity and is almost constant. 


2.2 White dwarf cooling rate 

Now, we recapitulate the discussion of white dwarf cooling as outlined by Mestel ( 1952l 'l and Schwarzschild (195^). The 
thermal energy of the ions is the only significant source of energy that can be radiated when a star enters the white dwarf 
stage. As most of the electrons occupy the lower energy states in a degenerate gas, the thermal energy of the electron gas 
cannot be readily released. Also, the energy release from neutrino emission is considerable only in the very early phase when 
the temperature is high. The thermal energy of the ions and the rate at which it is transported to the surface and to be 
radiated are dependent on the specific heat, which in turn depends significantly on the physical state of the ions in the core 
as discussed below. 


2.2.1 Crystallization and melting temperature 

The discussion here is applicable to stars with ionic densities less than 10® g cm“® as is the case for the white dwarfs we study. 
At very high temperatures, T 10^ K, and the ions are in a gaseous state with the kinetic energy of thermal motion large 
compared to the Coulomb energy. The specific heat per ion c„ is thus that of an ideal gas with c„ = Sfcs/2. As the ionic core 
cools, the gas starts deviating from its ideal behavior and the Coulomb energy can no longer be neglected. In the absence of 
any thermal motion, the configuration of the ions that minimizes the Coulomb energy is a lattice (with a BCC structure). 
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The ions here are no longer free like in an ideal gas and their kinetic energy corresponds to vibrations about the equilibrium 
positions of the lattice. An estimate of the temperature Tg below which the ions can no longer be thought of as being in a 
gaseous state is obtained by equating the kinetic energy with the vibrational energy yielding 

Tg « 3 X K, (8) 

where Z is the atomic number. Below Tg, the ionic system can be thought of as being in a liquid state. As the core cools 
further, the liquid crystallizes releasing its latent heat of melting which temporarily arrests further cooling. However, once 
this energy has been radiated, the cooling starts again and the ions eventually crystallize at a temperature Tm, which is the 
melti ng temperature of the Coulomb lattice. An estimate of Tm can be obtained from the Lindemann criterion (ILindemann 
1910li . which states that melting occurs when the root mean square displacement of the ions about their equilibrium positions 


due to thermal vibrations is about a quarter of the inter-ionic spacing. This gives 




2 X 


■5/3 


K. 


( 9 ) 

Thermal vibrations of the ions correspond to excitations of the normal (phonon) modes of the ionic crystal. Below Tm, the 
distribution of energy in these modes can initially be thought of as being consistent with classical equipartition yielding 
= 3fcs, which is the Dulong-Petit law. Upon further cooling, the quantization of the energy in each mode needs to be 
considered, and the specific heat is that of a quantum crystal. The temperature below which this happens is the Debye 
temperature 9d oi the lattice and is related to the the plasma frequency flp of the ions as 

h 


Od = 


ks 


- Ht, 


where 


mi 


2„2 \ 1/2 


( 10 ) 


( 11 ) 


Here rn is the number density of the ions, m; is the mass of a single ion, e is the charge of an electron and h is Planck’s 
constant. For T 6n, 

4 / ^ X 3 




ks- 


( 12 ) 


2.2.2 Time dependence of temperature 


The c ooling rate of a white dwarf is given by —dU / dt, which can be equated to the expression for L to give (IShapiro fc Teukolskv 
1983ll 

d 


dt 


/ 


CvdT=CAmuT’^^^, 


where CMq « 2 x 10® erg s ^ and A is the atomic mass. For T ^ Tg, Cv ZkB/'i., which gives 


_d f ikBT/2 \ 

dt \ Amg J 


= ct'^''^. 


Integrating, we obtain 


3 ks 


5 Amu 


[T“5/2 _ ) = C{t - to) = Ct, 


(13) 


(14) 


(15) 


where Tq is the initial temperature (before cooling starts), T is the present temperature at time t and t = t — to is the age of 
the white dwarf. Using Eqn. (5) we rewrite Eqn. (15) as 


3 ksM 
5 Amu 


h 

Lo 


3 fcsM T 


5 Amg L 




(16) 


where, Lo and L are the initial and present luminosities of the white dwarf, respectively. Putting in the values of the constants 
in Eqn. (15), we have 

y-5/2 _ y-5/2 ^ 2.4058 x 10 ”®'‘t. (17) 

Using Eqn. (17), we find the present temperature, T, and luminosity, L, at the interface for different values of r which 
correspond to the present age of the white dwarf. We calculate T for T* given in Table 1 and r = 10 billion years, i.e., 
3.1536 X lO^’^ s. It is important to note that r cannot exceed 13.8 billion years, which is the present age of the universe. 

From Fig. 3, it can be seen that cooling at the interface is considerable only for higher values of luminosities (L ^ 10“®Lq). 
The value of T* remains almost constant for L — 10“®!/©. Using Eqn. (5) in Eqn. (16), we obtain 

'L-®/"-Lo“®/"' 


M-5/7 


(18) 
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Figure 3. Variation of temperature with time for different ini¬ 
tial luminosities: 10 “®-Lq (dashed line), 10 “^-Lq (dotted line), 
1 O“®-L 0 (dot-dashed line) and 10~^Lq (solid line). 



Figure 4. Variation of luminosity with time for different initial lu¬ 
minosities: 10 “®Lq (dashed line), 10~'^Lq (dotted line), 1 O“®L 0 
(dot-dashed line) and 1 O“^L 0 (solid line). 


Therefore, the luminosity and interface temperature both reduce with time, for a given initial luminosity of the white dwarf. 
From Fig. 4, it can be seen that even after 10 billion years, L decreases only by one order of magnitude. This explains why 
many white dwarfs have yet not faded from view, although their initial luminosities may have been quite low. Figure 3 shows 
that white dwarfs spend most of the time near their present temperature. 


For < T < T„, c„ 


3fcs, which gives, 

6 ksM 

5 Atria 


n 

Lq 


6 fcflM T 
5 Ania L 


1 - 


rp \ 5/2 


This r is bigger than that in Eqn. (16) by a factor of 2, because of the inclusion of the lattice potential energy. 
For T <C 0D, using Eqn. (12), we obtain 


1287r^ 


rp \ 1/2 

-1 


MfcsT 1287r^ ksM 


ArriaL 


Lo 


T 


(19) 


( 20 ) 


5 AniaOD 

This is shorter than the classical result whenever T ^ QA6d ( Mestel fc RudermanlllQOT l. In this case, c„ falls off rapidly below 
9d, which leads to faster cooling and better agreement with observations. 

We have neglected neutrino emission from the right-hand side of Eqn. (13 ). Thermal neutrino emission dominates over 
photon emission if L ^ 10~°'^Lq and T ^ 10^'® K I Shapiro fc Teukolskv 19831 ). However, the luminosities that we consider 
here satisfy L ^ 10~^Lq and hence, photon emission dominates over thermal neutrino emission. _ 


_ Moreo ver, convection might also resu lt in shorter cooling timescales due to more efficient energy transfer. But iLamb fc Van Horn 

1 1975h and Fontaine fc Van Horn 1 19761) showed that convection is not significant to a first approximation. 


3 TEMPERATURE PROFILE IN THE PRESENCE OF A MAGNETIC FIELD 


Here we solve the magnetostatic equilibrium and photon diffusion equations in the presence of a magnetic field (B) and 
investigate the change in the temperature profile inside a white dwarf. We perform our calculations for both radially constant 
and radially varia ble magnetic fields. Entire c alculation is being performed in the same spirit as done for non-magnetised 
white dwarfs (e.g. Shapiro fc Teukolskv 19831 ) above, except with the inclusion of magnetic fields. Hence, the hydrostatic 
equilibrium is replaced by the magnetostatic equilibrium. 

The presence of B inside a white dwarf gives rise to an additional magnetic pressure, Pb = which adds up with 

the electron degeneracy pressure to give rise to the total pressure. Furthermore, the density also has a contribution from 

□ 2 

the magnetic field which is given by, ps = Such a situation can be tackled more ingeniously in the general relativistic 
framework rather than Newtonian framework. Nevertheless, here as a first approximation, we construct the magnetostatic 
equilibrium and photon diffusion equations in a Newtonian framework approximately as 


d, 7-> \ GAd , , 

— (P + Pb) =- —(p + Pb), 

dr r‘‘ 


and 


dr 


3 ^0(p T Pb) L 


4ac 


ye. 5 


47rr2 ’ 


( 21 ) 


( 22 ) 
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Figure 5. Variation of density with temperature for B = 10^^ G 
and different luminosities: 1 O“^L 0 (dashed line), 10 “^^Lq (dotted 
line) and 10“^ Lq (dot-dashed line), p* and T* are obtained from 
the intersection of p — T profiles with Eqn. (4) (solid line). 



Figure 6. Variation of density with temperature for L = 10' -5Lq 
and different magnetic fields: 10^^ G (dashed line), 10^^ G (dotted 
line) and 2 x 10^^ G (dot-dashed line), p* and T* are obtained 
from the intersection of p — T profiles with Eqn. (4) (solid line). 


respectively, neglecting magnetic torsion terms. Using Eqns. (21) and (22), we can write 


d 


dT 




\ 4ac 4:7rGM 


Stt J 


kqL 


(^+ sf^) 


(23) 


It has already been shown that including th e effects of a magnetic pressur e gradie nt and magnetic density gives rise to stable 
highly super-Chandrasekhar white dwarfs ( Das &: Mukhopadhvav 2014a 20151 ). Note that a large number of magnetised 

white dwarfs with surface fields as high as 10® G have been discovered by the Sloan Digital Sky Survey (SDSS) ( Schmidt et al. 


20031) . It is possible that their central field is several orders of magnitude larger than the surface field. To capture the variation 


of B, we use a profile proposed bv iBandvopadhvav et al.l (119971 ). modeling B as a function of p, given by 

+ Bo 


B I ] ^B, 

mo. 


1 - exp -p — 
Po 


(24) 


where Bs is the surface magnetic field, Bq is a parameter having dimension of B, 77 and 7 are parameters determining how 
the magnetic field decays from the core to the surface. The value of po is chosen to be ten percent of pc. We set, 77 = 0.8, 
7 = 0.9 and po = 10® g cm“® for all our calculations. Close to the surface, we have p —>■ 0 and therefore B ^ Bg. 

In §3.1, we first find the temperature profile in the presence of a constant radial magnetic field by setting Bo = 0 and 
choosing various values of Bs, so that the magnetic field is independent of p. Then in §3.2, we vary both Bo and Bs to find 
the temperature profile for a radially varying field. 


3.1 Constant magnetic field 


By constant magnetic field, we mean a magnetic field independent of p and having the same magnitude throughout the white 
dwarf. This is somewhat of an ideal case, but is a good choice to start the investigation. For this, we set Bo — 0 which 
gives B — Bs = constant and choose L such that 10 ~®Bq ^ L 10~^Lq is satisfied. We consider a white dwarf with mass 
M = Mq and vary the magnitudes of Bs and L to check how T*, r« and the temperature profile in the envelope change. 

As we are interested in the surface layers which are nondegenerate, we can substitute for P in terms of p in Eqn. (23) by 
using the ideal gas equation of state to obtain 


d / pkBT\ _ 


4ac 47 rGM 


dT \ J 


k.oL 


{p+^) 


Putting in the values of the constants and using Z = 0.1 and A = 0, we have 

1.168 r® 


5.938 X lO'^p + 5.938 x 


(P 2 . 261 x 10 ^^ 


(25) 


(26) 


From Eqn. (22), we have 


^ = -8.758 X 10“®'^ 
dT 


{P+^) 


(27) 
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Figure 7. Variation of radius with temperature for B = 10^^ G 
and different luminosities: 10 “^Lq (dashed line), 10 -4 Lo (dotted 
line), 10“^ Lq (dot-dashed line) and 10”^ Lq (solid line). 



Figure 8. Variation of radius with temperature for L = 10“^ Lq 
and different magnetic fields: 10^^ G (dashed line), 10^^ G (dotted 
line), 5 x 10^^ G (dot-dashed line) and 10^^ G (solid line). 


Table 2: Interface temperature (in K) for different luminosities and magnetic fields 


L 

10 -® Lq 

O 

1 

o 

10"® Lq 

10"® Lq 

B = 0 

2.332 X 10® 

4.50238 X 10® 

8.69273 X 10® 

1.6783 X 10’’ 

B = 10“ G 

2.33202 X 10® 

4.5024 X 10® 

8.69275 X 10® 

1.6783 X 10’’ 

B = 10^^ G 

2.33478 X 10® 

4.50439 X 10® 

8.69418 X 10® 

1.67841 X 10^ 

B = 10“ G 

2.5559 X 10® 

4.68586 X 10® 

8.83246 X 10® 

1.68858 X 10’’ 

B = 5 X 10“ G 

3.92503 X 10® 

6.41893 X 10® 

1.07706 X lO’' 

1.8764 X lO’’ 

B = 10“ G 

5.07388 X 10® 

8.14314 X 10® 

1.32226 X lO’' 

2.19241 X lO’’ 

B = 5 X 10“ G 

9.54829 X 10® 

1.51536 X 10^ 

2.40817 X lO’' 

3.83716 X lO’’ 

B = 10“ G 

1.25884 X lO’’ 

1.99615 X 10^ 

3.16693 X lO’' 

5.02952 X lO’’ 

B = 5 X 10“ G 

2.39528 X lO’’ 

3.79647 X 10^ 

6.01765 X lO’' 

9.53938 X lO’’ 


We solve Eqns. (26) and (27) simultaneously by specifying boundary conditions at the surface as 

p(Ts) = Ps = 10"“ g cm■^ r{Ts) = R = 5000 km, Ts = 5000 K. (28) 

Once we obtain p — T and r — T profiles for the given boundary conditions, we can find T* and p* by solving p — T profile 
along with Eqn. (4), as shown in Figs. 5 and 6. This works because the p — T profile is valid in the whole envelope whereas 
Eqn. (4) is valid only at the interface. Once we know T*, we can also find r* from the r — T profile, as shown in Figs. 7 and 8. 

In Tables 2, 3 and 4, we present the values of T*, p* and r* respectively for different values of B and L. The values of 
T*, p* and r* have been specihed to high precision (up to higher decimal places), as the temperature after cooling is quite 
sensitive to the initial conditions. We see that T* and p* increase with both B and L. For a given B, T* becomes almost twice 
in magnitude as L increases by one order of magnitude. Also, r* decreases with an increase in L for a given B. Therefore, the 
interface moves inwards with increase in L for a given B. 

Table 3: Interface density (ingcm"^) for different luminosities and magnetic fields 


L 

10"® Lq 

0 

1 

o 

10 ® Lq 

10 ^ Lq 

B = 0 

170.722 

457.995 

1228.66 

3296.13 

B = 10“ G 

170.725 

457.998 

1228.67 

3296.13 

B = 10“ G 

171.028 

458.302 

1228.97 

3296.43 

B = 10“ G 

195.891 

486.275 

1258.41 

3326.46 

B = 5 X 10“ G 

372.789 

779.635 

1694.55 

3896.59 

B = 10“ G 

547.908 

1114 

2305.01 

4921.32 

B = 5 X 10“ G 

1414.45 

2827.94 

5665.37 

11395 

B = 10“ G 

2141.17 

4275.51 

8543.9 

17099.7 

B = 5 X 10“ G 

5619.96 

11214.2 

22378.9 

44666.1 
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Luminosity and cooling of highly magnetised white dwarfs 9 


As L increases, dT/dr near the surface increases. The value of dT/dr near the surface also increases with the increase in 
B whereas r — T relation is almost linear, as for the non-magnetised case (or for B ^ 10® G). 


Table 4: Interface radius (in R) for different luminosities and magnetic fields 


L 

10 -® Lq 

10-^ Lq 

10 ® Lq 

10-® Lq 

B = Q 

0.9978 

0.9958 

0.9918 

0.9844 

B = 10^° G 

0.9979 

0.9959 

0.9920 

0.9846 

B = 10^1 G 

0.9982 

0.9963 

0.9926 

0.9855 

B = 10^® G 

0.9991 

0.9978 

0.9949 

0.9888 

B = 5 X 10^® G 

0.99975 

0.999256 

0.99788 

0.994272 

B = IQi® G 

0.99987 

0.9996 

0.9988 

0.9965 

B = 5 X 10^® G 

0.999972 

0.999914 

0.999732 

0.999166 

B = lO^'^ G 

0.999986 

0.999958 

0.999866 

0.999576 

B = 5 X 10^“^ G 

0.999998 

0.999992 

0.999972 

0.999914 


3.2 Varying magnetic field 


It is generally believed that the magnetic field strength at the surface of a white dwarf is several orders of magnitude smaller 
than the central field strength. Therefore, in this section we consider a more realistic density dependent magnetic field 
profile such that the magnetic field strength decays from the core of the white dwarf to its surface. We choose L such that 
lO“®Z /0 ^ 1/ ^ 1 O“®L 0 and vary the magnitudes of Bs and Bq, keeping r/ and 7 constant, to investigate how T*, r« and 
temperature profile change. 

The calculations done in this section are very similar to those of §3.1 except that ^ j ^ 0 in Eqn. (23). As earlier, 
we substitute P in terms of p in Eqn. (23) using the ideal gas equation of state to obtain 

d f pkeT B‘^\ _ 4ac 47rGM 
dT \ Stt / 3 kqL 


Putting in the values of the constants and substituting Z = 0.1 and A = 0, we obtain 

5.938 X IOV +5.938 x + 0.07965^-^ = 

dT dp dT 


rjnb 


(yP 2 . 261 x 10 ^^ ) 


(30) 


Eqns. (27) and (30) are simultaneously solved by specifying boundary conditions at the surface as before: p{Tb) = 10“^°g cm“® 
and r{Ts) = R = 5000 km, where Ts = 5000 K. The values of T*, p* and r, are obtained from the p — T and r — T profiles 
in the same way as in the case of constant field, as shown in Figs. 9, 10, 11 and 12. Tables 5, 6 and 7 show T,, p* and r, 
respectively, for various values of 73 = (73s, Bq) and L. 


Table 5: Interface temperature (in K) for different luminosities and magnetic fields 


L 

10-® Lq 

10-" Lq 

10-® Lq 

10 ® Lq 

B = (0,0) 

2.332 X 10® 

4.50238 X 10® 

8.69273 X 10® 

1.6783 X 10 '^ 

B = (10® G, 10^^ G) 

2.44368 X 10® 

4.61565 X 10® 

8.81325 X 10® 

1.69275 X 10 '^ 

B = (10® G,5 X 10^"^ G) 

2.81701 X 10® 

5.1498 X 10® 

9.65885 X 10® 

1.84453 X 10^ 

B = (10^° G,10i® G) 

2.43875 X 10® 

4.60434 X 10® 

8.78775 X 10® 

1.68706 X 10'^ 

B = (10^° G,10^'‘ G) 

2.98466 X 10® 

5.25529 X 10® 

9.50633 X 10® 

1.76262 X 10 '^ 

B = (10^° G,5 X 10^^ G) 

3.93081 X 10® 

6.6721 X 10® 

1.15646 X 10 '^ 

2.06066 X 10 '^ 

B = (10“ G,10^® G) 

2.98321 X 10® 

5.24806 X 10® 

9.48594 X 10® 

1.75755 X 10 '^ 

B = (10“ G,10^'‘ G) 

4.4574 X 10® 

7.43835 X 10® 

1.25587 X 10 '^ 

2.15789 X 10 '^ 

B = (10“ G,5 X lO^'^ G) 

6.21249 X 10® 

1.02625 X 10^ 

1.70449 X 10 '^ 

2.85704 X 10^ 

B = (10^® G,10i^ G) 

7.3182 X 10® 

1.19236 X 10 '^ 

1.95736 X 10'^ 

3.23512 X 10 '^ 

B = (10^® G,5 X lO^'^ G) 

1.02447 X 10^ 

1.6724 X 10^ 

2.74177 X 10^ 

4.51058 X 10^ 

B = (10^^ G,5 X lO^'^ G) 

2.10308 X 10^ 

3.18665 X 10^ 

4.91886 X 10^ 

7.74337 X 10^ 


From Table 5, we see that T* increases with increasing 73s, Bq and L. For a given (73s, 73o), 7* increases as L increases. 
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Figure 9. Variation of density with temperature for B = 
(Bs^Bq) = (10^^ G, 10^^ G) and different luminosities: 10“^ Lq 
(dashed line), 10““^ Lq (dotted line) and 10“^ Lq (dot-dashed 
line), p* and T* are obtained from the intersection of the p — T 
profiles with Eqn. (4) (solid line). 



Figure 11. Variation of radius with temperature for B = 
(10^^ G, 10^^ G) and different luminosities: 10“^ Lq (dashed line), 
10“^^ Lq (dotted line), 10“^ Lq (dot-dashed line) and 10“^ Lq 
(solid line). 



Figure 10. Variation of density with temperature for L = 
10“^ Lq and different magnetic fields: (10^^ G, 10^^ G) (dashed 
line), (1012 G, G) (dotted line) and (10^2 G, 5 x lOi^ G) (dot- 
dashed line), p* and T* are obtained from the intersection of the 
p — T profiles with Eqn. (4) (solid line). 



Figure 12. Variation of radius with temperature for L = 10~^ Lq 
and different magnetic fields: (IQH G, IQI'I G) (dashed line), 
(10i2 G, G) (dotted line), (IO 12 G, 5 x IQi'i G) (dot-dashed 
line) and (IQl^ G, 5 x G) (solid line). 


However, the percentage change in T* with the change in L decreases as Bs and Bq increase. In other words, the increase of 
T* due to the increase of L is somewhat saturated with the increase in 5 = (5s, Bq). Also, for a given Bs and 5o, r* decreases 
with L. Therefore, the interface moves inwards with an increase in L for a given (5s, 5o). 


Table 6: Interface density (in g cm ^) for different luminosities and magnetic fields 


L 

10-® Lq 

10-^ Lq 

10-® Lq 

10-" Lq 

B = (0,0) 

170.722 

457.995 

1228.66 

3296.13 

B = (10® G,10^^ G) 

183.132 

475.387 

1254.3 

3338.76 

B = (10® G,5 X 10^"^ G) 

226.663 

560.252 

1439.08 

3797.74 

B = (10^° G, 10^® G) 

182.578 

473.641 

1248.86 

3321.95 

B = (10^° G, 10^'* G) 

247.196 

577.555 

1405.13 

3547.6 

B = (lO^® G,5 X lOi'* G) 

373.612 

826.213 

1885.37 

4484.42 

B = (10“ G, 10^® G) 

247.015 

576.363 

1400.61 

3532.32 

B = (10“ G, 10^'‘ G) 

451.149 

972.552 

2133.61 

4805.54 

B = (10“ G,5 X lO®'^ G) 

742.329 

1576.09 

3373.56 

7321.02 

B = (10®^ G, 10®'‘ G) 

949.083 

1973.82 

4151.48 

8821.32 

B = (10^® G,5 X 10^^ G) 

1571.98 

3278.76 

6882.48 

14522.7 

B = (10^® G,5 X 10^^ G) 

4623.6 

8623.83 

16538.4 

32665.8 
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B = (10^2 G, lO^'* G) 

0.7492 

0.7158 

0.6851 

0.6560 

B = (10^2 G 5 X 1014 G) 

0.3485 

0.3165 

0.2907 

0.2694 

B = (10^® G,5 X lO^"* G) 

0.2512 

0.1892 

0.1454 

0.1152 


As shown in Figs. 11 and 12, unlike the constant magnetic field case the r — T profile is no longer linear for any L. Also, 
as L increases, dT/dr near the surface increases. The value of dT/dr also increases more with the increase in Bs as compared 
to increase in Bq. dT/dr near the surface is very high for Bs ^ 10^^ G and Bo ^ 10^® G. Therefore, the temperature fall near 
the surface increases for increasing strength of the magnetic field and luminosity, as for the previous cases. The values of p* 
and r* also increase and decrease respectively with the increase of L and/or B. 


4 VARIATION OF LUMINOSITY WITH MAGNETIC FIELD 

In this section, we determine how the luminosity of a white dwarf changes as the magnetic field strength increases such that, 

1 . the interface radius for a magnetised white dwarf is the same as that for a non-magnetised white dwarf i.e., = T'*,fl=o, 

2 . the interface temperature for a magnetised white dwarf is the same as that for a non-magnetised white dwarf i.e., = 

Tt,B=0- 

We also check if the conditions above can be satisfied simultaneously by changing the value of Ts accordingly. We calculate 
L for various magnetic field profiles, such that either r* or T* has the same value as that for the non-magnetised white dwarf 
with a luminosity L — 10~®Lq. 


4.1 Fixed interface radius 
4-1.1 Constant magnetic field 

Note that for B = 0 and L = 10~® Lq, we have found r* = 0.9978 R, p* = 170.722 gcm“® and T* = 2.332 x 10® K (from 
Table 1). We solve Eqns. (26) and (27) as done in §3.1, but this time we vary L in order to fix r* = 0.9978 R. We use the 
same boundary conditions as in Eqn. (28). 

Table 8: Variation of luminosity with magnetic field for a fixed r* 


B (in G) 

L (in Lq) 

T* (in K) 

p. (in g cm ®) 

B = 0 

10 "® 

2.332 X 10® 

170.722 

B = 10’^ 

1.045 X 10“® 

2.36151 X 10® 

174.401 

B = 10^® 

1.173 X 10”® 

2.4408 X 10® 

182.808 

B = 10“ 

1.882 X 10”® 

2.79632 X 10® 

224.17 

B = 10^2 

1.05 X lO""' 

4.74811 X 10® 

495.998 

B = 5x 10^^ 

1.092 X 10”® 

1.09932 X 10'^ 

1747.37 

B = 10^® 

3.616 X 10"® 

1.74767 X 10'^ 

3502.58 

B = 5x 10^® 

7.305 X 10”^ 

5.76089 X 10'^ 

17473.7 


Table 8 shows that the luminosity has to increase with an increase in magnetic field strength if the interface is at the same 
radius. However, the increase in L is appreciable only for B > 10^^ G. The value of T* also increases with an increase in B, 
which is expected as T* is proportional to L. 


4-1-2 Varying magnetic field 

Here, we assume a magnetic field profile as given by Eqn. (24) and find the variation of luminosity with a change in Bs and 
Bo so that the interface radius is same as for the non-magnetic case. 

We solve Eqns. (27) and (30) as done in §3.2, but this time vary L in order to fix r* = 0.9978 R. 

Table 9: Variation of luminosity with magnetic field for a fixed r* 


B (in G) 

L (inL©) 

T* (in K) 

p* (in gem ®) 
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B = (0,0) 

10 “® 

2.332 X 10® 

170.722 

B = (10®, 10^^) 

9.78 X 10"® 

2.3184 X 10® 

170.722 

B = (10®, 10^®) 

5.17 X 10"® 

1.94342 X 10® 

129.882 

B = (10®, 5 X 10^®) 

2.87 X 10"® 

494967 

16.6941 

B = (10“,10^®) 

8.15 X 10"® 

2.21142 X 10® 

157.654 

B = (10“,5 X 10^®) 

1.35 X 10"® 

1.37434 X 10® 

77.2393 

B = (10“,10i®) 

3.08 X 10"® 

542023 

19.1304 

B = (lO^SlO^®) 

2.33 X 10"®® 

70709.3 

0.901391 


Interestingly, Table 9 shows that L and T* both decrease as the magnetic field strength increases. However, the change is 
appreciable only for Bg ^ 10^° G and Bo ^ 10^^ G. The value of L becomes quite low ~ 10“® Lq and lower, for white dwarfs 
with (Bs, Bo) = (10^° G, 5 X 10^^ G) and higher, which might make it difficult to detect such highly magnetised white dwarfs. 


4.2 Fixed interface temperature 

4-2.1 Constant magnetic field 

First we look at white dwarfs with constant a magnetic field strength throughout. We solve Eqns. (26) and (27) as done in 
§3.1, but this time we vary L in order to fix T* = 2.332 x 10® K. We use the same boundary conditions as in Eqn. (28). 

Table 10: Variation of luminosity with magnetic field for fixed T* 


B (in G) 

L (in Lq) 

p* (in g cm ®) 

r* 

B = 0 

10 ’® 

170.722 

0.9978 R 

B = 10® 

10 ”® 

170.722 

0.9978 R 

B = 10“ 

9.9583 X 10"® 

170.723 

0.9982 R 

B = 10“ 

6.9672 X 10“® 

170.723 

0.9992 R 

B = 5 X 10“ 

8.079 X lO"’® 

170.722 

0.999926 R 

B = 10“ 

2.143 X 10"'^ 

170.722 

0.99998 R 

B = 5 X 10“ 

8.7423 X 10"® 

170.723 

R 

B = 10“ 

2.1869 X 10"® 

170.723 

R 


Table 10 shows that unlike the case with fixed r* and constant magnetic field, L here has to decrease with the increase in 
B. However, the decrease in L is significant only when B ^ 5 x 10^^ G. For higher magnetic field strengths, the luminosity 
becomes very small and goes beyond the detection limit. Also, r* gradually increases as B increases. This is because the rate 
of temperature fall with radius across the envelope is higher for higher fields (from Eqn. 27). Therefore, r* has to increase for 
T* to be constant. 


4-2.2 Varying magnetic field 

Here, we solve Eqns. (27) and (30) as done in §3.2, but this time we vary L to get T* = 2.332 x 10® K, using the same 
boundary conditions as in Eqn. (28). 

Table 11: Variation of luminosity with magnetic field for a fixed T* 


B (in G) 

L (in Lq) 

p* (in g cm ®) 

r* 

B = (0,0) 

10 ~® 

170.722 

0.9978 R 

B = (10®,10“) 

9.9823 X 10"® 

170.722 

0.9978 R 

B = (10®,10“) 

8.4274 X 10”® 

170.722 

0.9940 R 

B = (10®, 5 X 10“) 

4.7561 X 10"® 

170.722 

0.97935 R 

B = (10“,10“) 

9.8251 X 10”® 

170.723 

0.9977 R 

B = (10“,10“) 

8.4848 X 10"® 

170.723 

0.9959 R 

B = (10“,10“) 

3.5514 X 10"® 

170.722 

0.97857 R 

B = (10“,5 X 10“) 

9.721 X 10"^ 

170.722 

0.90875 R 
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B = (10^\10^®) 

3.5439 X 10"® 

170.722 

0.9896 R 

B = (10^\10^'‘) 

5.1813 X 10"'^ 

170.723 

0.9192 R 

B = (10i\5 X 10^'*) 

1.07855 X 10"^ 

170.723 

0.69839 R 

B = 

3.5327 X 10"® 

170.723 

0.8414 R 

B = (10l^5 X 10i'‘) 

7.0921 X 10"® 

170.723 

0.51559 R 


Like the constant field case, we find that even for varying magnetic field strength, L has to decrease as B increases for T* to 
be unchanged. From Table 11, we see that L becomes very small when Bs > 10^° G and Bo ^ lO^'^ G. We also see that r* 
decreases with increase in magnetic field strength. 


4.3 Fixed interface radius, fixed interface temperature and varying surface temperature 

We find that the p — T and r — T profiles are sensitive to Ts only for relatively low values of the magnetic field: 

1. B ^ 10^° G for constant magnetic field. 

2. Bs ^ 10^° G and Bq ^ 10^^ G for a varying magnetic field. 

Hence, for the larger values of B that we are considering, we do not have any considerable change in the values of T* and r* 
(as obtained for Ts = 5000 K), by changing Ts. 


5 COOLING IN THE PRESENCE OF A MAGNETIC FIELD AND POST COOLING TEMPERATURE 
PROFILE 

We know how to find the cooling timescale of a white dwarf if we know the L — T relation, which was done in §2.2. In this 
section, we find the L — T relations for the cases mentioned in §4 and using those relations, implement cooling with time 
to find the present interface temperature, T*^pr from the initial interface temperature, r*,in, for r = 10 billion years. Once 
we know Tt,pr, we can easily calculate the present surface temperature, Ts,pr, using the equation for dT/dr relation and the 
present value of the luminosity Z/*,pr. 


5.1 Specific heat in the presence of magnetic field 


A magnetic field can, in principle, affect the state of the ionic core and thus its thermodynamic properties, such as the specific 
heat. The relevant parameter to quantify this effect is 

(31) 


b — 

(jJr) 


where 


ZeB 


ATvZ'^e^n 


are the ion cyclotron and ion plasma frequencies, respectively. Here, n is the number density of the ions. As discussed in §2.2, 
Up is effectively the Debye frequency of the ionic lattice. We would expect the effect of the magnetic field on the ionic core to 
be strong when 6 1, i.e. when the cyclotron frequency is comparable to or large r than t he De bye frequency of the lattice. 

The effect of a magnetic field on a BCG Coulomb lattice has been studied by Baikol ( 20091 ) . who concluded that there is 
an appreciable change of the specific heat only for 6 1 except when T Od- For the white dwarfs we consider, B < 10^^ 

G at the interface, which corresponds to fc Si 1. Further, the interface temperature is not significantly smaller than 6d. Thus, 
we are justified in working with a value of the specific heat appropriate for a non-magnetised system despite the presence of 
a magnetic field. 

It will be of interest to study the effect of much stronger magnetic fields on the ionic core and its specific heat. In 
particular, if the magnetic field is strong enough to cause Landau quantization of the electron gas in the core, it could result 
In a change in the effective ion-ion interaction as mediated by the electrons. This would be in addition to the direct effect of 
the field on the ionic core as described above. The effect of a magnetic field on the phonon spectrum of ions in conventional 
solid state systems has been investigated and found to be weak for practical values of the field appropriate to these systems 
I Holz 1972lf . However, in an astrophysical context, the effect might be appreciable if fields of the order of 10^® G arise and 
could result in very interesting physics. We defer an investigation of this to a future work. 

Using the values of T* for different L as obtained in §3, we can approximate L as a function of T*. Once we know how 
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L varies with T*, we can repeat the cooling calculations as done in §2.2 to first find (after cooling) and then hnd Ts,pr, 
using the value of in the dT/dr relation. 


5.2 Fixed interface radius 

5.2.1 Constant magnetic field 

From Table 8, we know the values of Z/*,in and the L — T relation can be found using Table 2. The calculated values of 
and Ts.pr for r* = 0.9978 R are given in Table 12. 

From Table 12, we see that as B increases, the coefficient in the L — T relation decreases whereas the exponent of T 
increases. Also the cooling is enhanced for higher Bs. We see an appreciable reduction in T* for B > 10^^ G. Hence, the 
change in B affects T* considerably. 


Table 12: Change in T* with time due to the presence of a magnetic field for fixed r* and Ts = 5000 K 


B (in G) 

Tg.,in (in K) 

L*,i„ (in Lq) 

L (in terms of T) 

Tt,pr (in K) 

0 

2.332 X 10® 

10 “® 

2.013 X io®r® ®°® 

2.2232 X 10® 

10 ® 

2.36151 X 10® 

1.045 X 10'® 

2.013 X io®r® ®®® 

2.24814 X 10® 

10 ^° 

2.4408 X 10® 

1.173 X 10"® 

2.013 X io®r® ®®® 

2.31443 X 10® 

10 “ 

2.79632 X 10® 

1.882 X 10"® 

1.994 X io®r® ®°® 

2.59905 X 10® 

10 “ 

4.74811 X 10® 

1.05 X 10"'‘ 

7.911 X io®r® ®®® 

3.86949 X 10® 

5 X 10“ 

1.09932 X 10^ 

1.092 X 10“® 

1.467 X 10 

7.18397 X 10® 

10 “ 

1.74767 X 10’’ 

3.616 X 10”® 

1.355 X io-2t^-558 

1.01002 X 10’’ 

5 X 10“ 

5.76089 X 10’’ 

7.305 X 10"® 

1.265 X io-®T'* ®“ 

2.27239 X 10^ 


5.2.2 Varying magnetic field 

We find the L — T relations for different values of B using Table 5. We also know the values of L*,i„ for different field strengths 
from Table 9. Using these, we calculate T,^pr and then Ta,pr for different values of B, as given in Table 13. 

Table 13: The change in T* with time due to the presence of a magnetic field for a fixed r* and Tg = 5000 K 


B (in G) 

Tt,in (in K) 

L^^in (in Lq ) 

L (in terms of T) 

T*,pr(inK) 

(0,0) 

2.332 X 10® 

10 "® 

2.013 X io®r® ®®° 

2.2232 X 10® 

(10®,10“) 

2.3184 X 10® 

9.78 X 10"® 

2.002 X 10®T® ®®® 

2.21225 X 10® 

(10®,10“) 

1.94342 X 10® 

5.17 X 10"® 

1.905 X io®r® ®®® 

1.88494 X 10® 

(10®, 5 X 10“) 

4.94967 X 10® 

2.87 X 10"® 

2.622 X io®r® ®®® 

4.94599 X 10® 

(10“,10“) 

2.21142 X 10® 

8.15 X 10"® 

1.905 X io®r® ®®® 

2.12125 X 10® 

(10“,5 X 10“) 

1.37434 X 10® 

1.35 X 10"® 

1.533 X io®T® ®“ 

1.35709 X 10® 

(10“,10“) 

5.42023 X 10® 

3.08 X 10"® 

5.588 X io'‘t® “® 

541630 

(10“,10“) 

7.07093 X 10'‘ 

2.33 X 10"“ 

3.411 X lO^T^ ^® 

70709.3 


Unlike the constant field case, here L decreases with increasing B. With the increase of field strength, the coefficient in the 
L — T relation decreases whereas the exponent increases. Unlike the constant held case, increasing B results in slower cooling. 
However, as before, Tg is not affected much by the increase of held strength. 


5.3 Fixed interface temperature 
5.3.1 Constant magnetic field 

We hnd the L — T relations for different values of B using Table 2 and then use Table 10 to obtain the values of L^^in for 
different helds. The values of r,.pr and Ts.pr calculated for different helds are given in Table 14. 
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Table 14: The change in T* with time due to the presence of a magnetic field for a fixed T* and Ts = 5000 K 


B (in G) 

(inK) 

(in Lq) 

L (in terms of T) 

r*,pr. (in K) 

0 

2.332 X 10® 

10 ’® 

2.013 X io®r® ®®° 

2.223 X 10® 

10 '* 

2.332 X 10® 

10 ~® 

2.013 X io®r® ®®° 

2.223 X 10® 

10 “ 

2.332 X 10® 

9.9583 X 10"® 

1.890 X io®r® ®®'‘ 

2.224 X 10® 

10 “ 

2.332 X 10® 

6.9672 X 10"® 

1.502 X 

2.262 X 10® 

5 X 10“ 

2.332 X 10® 

8.079 X 10"'^ 

4.875 X 10"® r'‘ ®®® 

2.325 X 10® 

10 “ 

2.332 X 10® 

2.143 X 10"^ 

1.066 X 

2.330 X 10® 

5 X 10“ 

2.332 X 10® 

8.7423 X 10"® 

6.210 X i0"^r^ ®®® 

2.332 X 10® 

10 “ 

2.332 X 10® 

2.1869 X 10"® 

1.350 X 

2.332 X 10® 


From Table 14, we find that the coefficient in the L — T relation decreases and the exponent of T increases with the increase 
of B. Also the cooling rate decreases with the increase of field strength and T* is almost unchanged even after r = 10 billion 
years for B ^ 10^^ G. 


5.3.2 Varying magnetic field 

The L — T relations for different values of B are obtained using Table 5 and the values of L*,in for different fields are obtained 
from Table 11. We then calculate r*,pr and Ts^pr for the different values of B given in Table 15. 

As for constant field, an increase of the magnetic field strength results in a decrease in the coefficient and increase in the 
exponent of T in the L — T relation, as shown in Table 15. In this case as well, the cooling rate decreases appreciably with an 
increase in magnetic field strength for Bs 10“ G and Bo 10“ G. 

Table 15: The change in T* with time due to the presence of magnetic field for a fixed T* and Ts = 5000 K 


B (in G) 

r*,i„ (in K) 

L*,in (in Lq ) 

L (in terms of T) 

Tt^pr (in K) 

(0,0) 

2.332 X 10® 

10 "® 

2.013 X io®r® ®°® 

2.223 X 10® 

(10®,10“) 

2.332 X 10® 

9.9823 X 10"® 

2.002 X io®r® ®°° 

2.224 X 10® 

(10®,10“) 

2.332 X 10® 

8.4274 X 10"® 

7.102 X io®r® ®®® 

2.238 X 10® 

(10®, 5 X 10“) 

2.332 X 10® 

4.7561 X 10"® 

9.652 X lO^T®-®®® 

2.274 X 10® 

(10“, 10“) 

2.332 X 10® 

9.8251 X 10"® 

1.794 X io®r® ®°^ 

2.225 X 10® 

(10“, 10“) 

2.332 X 10® 

8.4848 X 10"® 

6.877 X io®r® ®®® 

2.237 X 10® 

(10“,10“) 

2.332 X 10® 

3.5514 X 10"® 

2.040 X io®r® ®“ 

2.291 X 10® 

(10“,5 X 10“) 

2.332 X 10® 

9.721 X 10"^ 

7.268 X i0"ir‘‘ ®®® 

2.323 X 10® 

(10“, 10“) 

2.332 X 10® 

3.5439 X 10"® 

1.309 X 

2.298 X 10® 

(10“,10“) 

2.332 X 10® 

5.1813 X 10"'^ 

1.753 X I0"2r'* ®®'‘ 

2.328 X 10® 

(10“,5 X 10“) 

2.332 X 10® 

1.07855 X 10"^ 

1.675 X i0"®r^ ®“ 

2.331 X 10® 

(10“,10“) 

2.332 X 10® 

3.5327 X 10"® 

1.588 X io-5t'^ s64 

2.332 X 10® 

(10“,5 X 10“) 

2.332 X 10® 

7.0921 X 10"® 

1.534 X lo-fiT"* ®!® 

2.332 X 10® 


5.4 Fixed interface radius, fixed interface temperature and varying surface temperature 

As discussed in §4.3, the p — T and r — T profiles are sensitive to Ts only at relatively low values of B. Hence, for larger 
values of B, we do not have any considerable change in the values of T* and r* (as obtained for Ts = 5000 K), by changing 
Ts. Therefore, we do not calculate the cooling rates for this case as the change in T* is observable only for smaller values of 
B, which is the same as the non-magnetic result. 


6 DISCUSSION 

In this section, we discuss our results described in the previous sections and their physical significance. 
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6.1 Non-magnetised white dwarfs 


From Table 1, we see that as L increases in the envelope, T* and p* both increase whereas r* decreases. Using Eqns. (4), (5) 
and ( 6 ), we obtain 


L oc T* ®, p* oc and (T* 




(33) 


Hence, a white dwarf having a larger T* has more stored thermal energy, which it can radiate thus giving rise to a larger L. 
Also, a larger T* corresponds to a larger p*. For a fixed Ts and R, r* should decrease as T* increases. This is because the 
outer regions of the white dwarf are cooler as compared to the inner ones. 

We also hnd that | AT/Ar|, i.e., KTs — r*)/(i? — r*)| and the cooling rate AT/At increase with the increase in luminosity 
of the white dwarf. From Eqns. ( 6 ) and (15), 


AT 

Ar 


f n-Ts \ 
\R-r, J 


and (T"®/^ - T“®/^) oc r. 


(34) 


Note that L corresponds to the energy flux which is transported across a spherical surface and hence a larger luminosity 
means a larger flux (for a given radius) i.e. a larger AT/Ar. From Eqn. (34) it is clear that hotter or more luminous white 
dwarfs cool faster as To is larger. Therefore, the cooling rate should be more for a white dwarf having larger luminosity. 


6.2 Magnetised white dwarfs having fixed interface radius 

In §4.1, we have found out how much the luminosity has to increase or decrease for a magnetised white dwarf for it to have 
the same interface radius as a non-magnetised white dwarf. Then in §5.1, we have also computed the cooling rates for the 
corresponding cases and used the values of L as obtained in §4.1 to obtain the surface temperatures. Here we discuss our 
results. 


6.2.1 Constant magnetic field 


For a hxed r, and constant B, we have computed T,^in, p*, T*,pr. and Ts,pr. We have fixed r* = 0.9978 R, corresponding to 
L = 10-® Lq for the non-magnetic case. We have found that as B increases, L, T* and p* all have to increase for r* to be 
fixed, which can be seen from Table 8 . From Eqn. (26), we have 


dp 1 

dr'^T 


1.168 T®-® 

^ {P+-^) 


5.938 X 10 



(35) 


From Eqn. (35), we see that if B increases, then dp/dT decreases for a given T and L. As the initial conditions are 
the same, we obtain a smaller value of p for a given T, if B is larger. Hence, the presence of magnetic held suppresses the 
matter density at a given temperature. As a result, T,, which is obtained by the intersection of the p — T curve with Eqn. 
(4), increases. This is because the density of matter at which it becomes degenerate is reached only at a higher temperature. 
Similarly, an increase of L also leads to the increase of T*. Therefore, the increase of B leads to the increase of T* and L. 
Also, p* increases as p* oc T*^^ by Eqn. (4). 

We hnd that the temperature falls with radial distance from the centre of the white dwarf. | AT/Ar |, and the rate of 
cooling of the white dwarf, i.e., AT/At, both increase with B. As T* increases with increasing B and L and r* remain hxed, 
an increase in | AT/Ar | is expected. We have L to be of the form aT^, where a and 7 are constants which depend on the 
magnetic held strength, as given in Table 12. With the increase of B, a decreases and 7 increases. Hence, we have 


_d / 3fcBT \ 
dt \ Amp ) 


= aT^, 


(36) 


(37) 


which gives 

_ 3fcB {T^-^ -TI-^) 

Amp 0(7 - 1 ) 

When B increases, 0(7 — 1) decreases but (T^ ' — T^ ’') decreases more such that their ratio decreases. This is because To 
increases more in comparison to T with increasing B, as can be seen from Table 12. Therefore, the cooling rate increases with 
the increase of B. 


6.2.2 Varying magnetic field 

In §4.1.2 and §5.2.2, we have hxed r* for a varying magnetic held and calculated T,pn, Tt,pr, p* and Ts,pr. We have used 
r* = 0.9978 R, i.e. the value of interface radius corresponding to H = 0 and L — 10~® Tq. From Table 9, we have seen that 


© 2015 RAS, MNRAS OOO.ITII201 






Luminosity and cooling of highly magnetised white dwarfs 17 


as B increases, L, T, and p, all decrease for r« to be fixed. This is the opposite of what is observed for a constant magnetic 
field. 

From Eqn. (30) we have, 


dp 

dT 




6.5 


- 5.938 X lO^p 


(^5.938 X 107T + 0.07965^^ 


(38) 


It can be seen that Eqn. (38) is quite similar to Eqn. (35), with the denominator having an additional term In Eqn. (38). For 
the B configuration that we have considered, the strength of the field increases with density. Therefore, BdB/dp is positive 
and we obtain a smaller value of dp/dT for a given field strength for the varying magnetic field case as compared to the one 
with a constant magnetic field. Therefore, a varying magnetic field suppresses the matter density at a given temperature more 
than a constant field and we obtain a larger value of T,, for a similar reason as the constant field case. However, it should be 
noted that r* also decreases such that AT/Ar decreases. From Eqn. (27), 

, 2 


dT 


{p+£^) 


dr re s r 2 ’ 

In this case, the decrease in p with the increase in B is large enough to counteract an increase in ps and dT/dr decreases 
with the increase in B for sufficiently small L. This does not happen in the constant field case as the reduction of p is not 
large enough and dT/dr near the interface increases with increasing B. Therefore, in this case we have a smaller T* and a 
smaller L for larger field strengths, for r, to be constant. The value of p* also decreases with field as p, oc 

We find that \AT/Ar\ and AT/At both decrease with B. As T* decreases with the increase in B and r* remains fixed, 
a decrease in |AT/Ar| is expected. As for the previous case, L is of the form aT'^ as given in Table 13. Hence, we have 

(yl-T _ jU-T) 


r oc 


(40) 


a(7 - 1) 

When B increases, 0(7 — 1) and (T^ — Tg”^) both decrease. However, the decrease in 0(7 — 1) is more such that r increases. 

With increasing B, Tq and 7 do not change considerably whereas a decreases by orders of magnitude. Therefore, the cooling 
rate decreases with the increase in B. 


6.3 Magnetised white dwarfs with fixed interface temperatnre 

In §4.2, we have computed the change of L of a magnetised white dwarf for it to have the same T* as a non-magnetised white 
dwarf. Then in §5.3, we have found the cooling rates for the corresponding cases and used the values of L as obtained in §4.2 
to obtain Ts- Here we discuss our results. 


6.3.1 Constant magnetic field 


For a fixed T* and constant magnetic field, we have computed Tt,,in, r*, p* and Tt,,pr. We have fixed T* = 2.332 x 10® K, which 
is the interface temperature corresponding to L — 10~® Lq for the non-magnetic case. We have found that as B increases, L 
decreases and r* increases whereas p* is constant, as can be seen from Table 10. From Eqns. (26) and (27), we have 


and 


dp 1 1 1.168 T®-® 

(p+s&) 


5.938 X 10 



dr 

dT 


rjnb 


{P+^) 


(41) 


(42) 


As p* oc for a nondegenerate envelope, p* has to be fixed because T* is fixed. We know from § 6.2 that the presence of a 
magnetic field suppresses the matter density for a given T. As the initial conditions for the p — T profile are the same, dp/dT 
has to be larger near the interface for the magnetic case, which is achieved by the reduction of L. Therefore, for T to be fixed 
with increasing field, L should decrease. 

Again as the initial conditions for solving Eqn. (42) are the same as for the non-magnetic case, dr/dT should be higher 
near the interface for obtaining higher r, for a given T*. This is achieved as the decrease in L is sufficiently small in comparison 
to the increase in pB such that dr/dT near the interface increases and we obtain a larger value for r*. An increase in |AT/Ar| 
is expected as T* is fixed whereas r« increases. 

We find that the cooling rate AT/At decreases with the increase of the magnetic field strength. As in §6.2, L can be 
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written as aT'^ as given in Table 14 and the expression for the cooling timescale is given by Eqn. (37). Unlike the fixed r* and 
constant field case, here the decrease in — Tq is not as much as that for 01(7 — 1) which makes r larger for larger B. 

As mentioned earlier in §6.2, a changes by orders of magnitude whereas Tq and 7 do not change considerably. 


6.3.2 Varying magnetic field 

In §4.2.2 and §5.3.2, we have fixed T* for a varying magnetic field and calculated r*, p* and Tt,pr. We have fixed 

T* = 2.332 X 10® K, which is the interface temperature corresponding to L = 10“® Lq for the non-magnetic case and found 
that as B increases, both L and r* decrease, whereas p* remains constant, as can be seen from Table 11. We recall Eqns. (38) 
and (39) which are applicable here as well. 

As p* cx for a nondegenerate envelope, p* has to be fixed as T* is fixed. Also, we know from § 6.2 that a magnetic 
field suppresses p for a given T. As the initial conditions for the p — T profile are same, we should have larger dp/dT near 
the interface for the magnetic case, which happens because of a reduction in L. Therefore, for T* to be fixed with increasing 
field, L has to decrease. 

Now, as the initial conditions for solving Eqn. (39) are the same as those for the non-magnetic case and dT/dr near the 
surface is larger for larger magnetic fields (from Fig. 12), dT/dr should be smaller near the interface for obtaining a smaller 
value of r* for a given T* in this case. We find that with increasing B, the luminosity for the varying field is sufficiently small 
in addition to p being small, which counteracts the increase in ps thus making dT/dr near the interface smaller. Therefore, 
r* decreases with increasing B for fixed T*. AT/Ar decreases as T* is fixed whereas r* decreases. 

We find that the cooling rate AT/At decreases as magnetic field strength increases. We have, L — aT'^ as given in Table 
15 and the expression for the cooling timescale is given by Eqn. (37). As for the constant field and fixed T* case, here the 
decrease in 0(7 — 1) is larger as compared to the decrease in (T^~^ — Tq *'), which makes r larger for larger B. 


6.4 Highly magnetised white dwarfs and their invisibility 


For the solar neighbourhood of radius Ro ~ 8.0 kpc and disk thickness t ~ 0.5 kpc, the expected number of visible white 
dwarfs with number density 1.5 x 10“^ pc“® ( Shapiro fc Teukolskvlll983l l is 

Nvis.wd = TT X i?o X t X 1.5 X 10”^ pc~® « 1.5 X 10®. (43) 

As the average mass of observed white dwarfs is 0.6 Mq, the mass of these white dwarfs in the solar neighbourhood is 

= 0.6 Mq X 1.5 X 10® = 9 X 10® Mq. (44) 

Now, 97% of the main sequence stars, with density in the solar neighbourhood 0.064M© pc~®, ultimately turn out to be white 
dwarfs. Hence the upper limit of total mass due to white dwarfs in the solar neighbourhood is given by 

Mtotai.wd — TV X Rq X t X 0.97 X 0.064 M© pc~® « 6.24 x 10® Mq. (45) 

From Eqns. (44) and (45), we have an approximate estimate of mass corresponding to non-visible white dwarfs 


^non — vis,wd — ^total.wd ^vis.wd — 5.34 X 10 AIq . 


(46) 


Interestingly, this is the same order as the missing mass, of density ~ 0.02 Mq pc in the solar neighbourhood, which is 


M„ 


= -K X Ro X t X 0.02 Mq pc ® « 2 X 10® Mq. 


Now, observed white dwarfs have luminosities in the range ~ 10“® Lq — 10~® Lq . In §4, we have shown that the luminosity 
of a magnetised white dwarf has to decrease for either r* or T* to be the same as a non-magnetised white dwarf, except for 
the case when there is a uniform magnetic field inside the white dwarf and r* is fixed. In other words, for a magnetised 
white dwarf to have similar interface parameters as a non-magnetised one, its luminosity needs to be much smaller - the 
luminosity can be as small as ~ 10“® Lq for very high magnetic fields. Indeed the surface temperatures of white dwarfs with 
different B are found/observe d neither to be coherently changing with B nor to be differing significantly for different B (e.g. 
Koester fc Chanmugam 199fll l. Such low luminosity white dwarfs are naturally not visible. Above estimate argues that such 
highly magnetised white dwarfs may be a significant component of Mnon-vis,-wd. Like neutron stars, they are supposed to be 
in the very bottom left of the H-R diagram and hardly visible. 


7 SUMMARY AND CONCLUSION 

We have investigated the effects of magnetic field in the luminosity and cooling of white dwarfs. To the best of our knowledge, 
this is the first systematic attempt towards this mission. This is very useful to account for observability of recently proposed 
highly magnetised white dwarfs, in particular the ones having central field 5 x G. However, we have deferred our 


© 2015 RAS, MNRAS 000.1111201 










Luminosity and cooling of highly magnetised white dwarfs 19 


investigation for white dwarfs with fields > 10^® G for a future work, which will affect the EoS significantly and, subsequently, 
might change the thermal conduction and, hence, observable properties more severely. 

We have computed the variation of luminosity of highly magnetised white dwarfs with magnetic field strength (both 
radially constant and varying cases) and evaluated the corresponding cooling timescales for white dwarfs having fixed interface 
radius or temperature. Except for the fixed interface radius with constant magnetic field cases, we have seen that the luminosity 
is suppressed with the increase in field strength, in addition to a reduction of cooling rates (marginally). Therefore, white 
dwarfs with higher magnetic fields have lower luminosities and slower cooling, for the interface radius or temperature to be 
fixed. 

For a similar gravitational energy (similar mass and radius), an increasing magnetic energy necessarily requires decreasing 
thermal energy for white dwarfs to be in equilibrium, which results in a decreasing luminosity. We have found that the 
luminosity could be as low as ~ 10“® Lq for a white dwarf with the central field ~ 5 x 10^^ G and the surface field 10^^ G, 
for a fixed interface temperature. As a result, such white dwarfs appear to be invisible with respect to current astronomical 
techniques. With decreasing the surface field, while the luminosity tends to be in the observable limit, it is still ~ 10“® Lq for 
the surface field 10® G. For a fixed interface radius, however, luminosity could be much lower, ~ Lq, for central and 

surface fields respectively ~ 10^^ G and 10^^ G. By decreasing surface field to ~ 10® G, it is still 10“® Lq, well below 
the observational limit. Therefore, such white dwarfs, while expected to be present in the universe, are virtually invisible to 
us, perhaps lying in the lower left corner in the H-R diagram. 
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